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Structure Functions and Parton Distributions 

Jianwei Qiu®^* 
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Ames, Iowa 50011, U.S.A. 


In this talk, I review the status of theoretical understanding of nuclear structure func¬ 
tions and parton distributions and discuss the constraints on nuclear parton distributions 
from existing data and the global QCD analysis. 

1. Introduction 

Much of the predictive content of perturbative QCD (pQCD) treatment of the hadronic 
hard scattering is contained in factorization theorems [ |l| . Their purpose is to separate 
long- from short-distance effects in scattering amplitudes. They supply perturbatively 
uncalculable long-distance effects with physics content in terms of well-dehned matrix 
elements, which allows them to be measured experimentally or by numerical simulation. 
They also dehne the normalization of short-distance factors, which allows them to be cal¬ 
culated perturbatively. Predictions follow when processes with different hard scatterings 
but the same nonperturbative matrix elements are compared. Thus, quark and gluon dis¬ 
tributions measured in deep inelastic scattering may be used to normalize the Drell-Yan 
or jet cross section. 

Deep inelastic lepton scattering has long been regarded as the cleanest probe of con¬ 
stituent substructure. Considerable interest, therefore, greeted the observation [ ^ by 
European Muon Collaboration (EMC) that the structure function F 2 {xb,Q‘^) of an iron 
nucleus differs in signihcant ways as a function of from that for deuterium. It was the 
discovery of the EMC effect that opened a door for systematic study of QCD dynamics 
in a nuclear environment, which has led to many new QCD phenomena, e.g., shadowing, 
saturation, and color glass condensate. 

In this talk, I review our abilities and limitations to generalize the pQCD factorization 
theorems to the hard scattering involving nuclei. Nuclear parton distributions are heavily 
used in phenomenological description of hard processes in heavy ion reactions at SPS 
and RHIC energies and in calculating predictions at the LHC energies. According to the 
factorization theorems, nuclear dependence of parton distributions should be universal 
or process independent. I discuss the constraints on nuclear parton distributions from 
existing data and the global QCD analysis. 


* jwqOiastate. edu. I thank Eskola and Salgado for useful discussions and figures. This work is supported 
in part by the United States Department of Energy under Grant No. DE-FG02-87ER40371. 



2 


2. Structure Functions and Parton Distribntions 
2.1. Strnctnre functions 

An inclusive deep inelastic scattering (DIS) process is generically of the form, i{E) + 
h{p) ^{E') + X, where £ represents a lepton, h a hadron (a nucleon or nucleus), and 
X an arbitrary hadronic hnal-state. The process, illustrated in Fig. 1(a), is initiated by 
the exchange of a virtual photon (or a vector boson in general). In DIS, the momentum 
transfer between the lepton and the hadron, q, is spacelike, —= Q^. The Bjorken 
scaling variable is dehned as = 2 ^^) where u is the energy transferred from the 

lepton to the hadron in the hadron (target) rest frame, v = E — E'. In the same frame, 
_ 4,EE' siv?{6/2) with the lepton scattering angle 6. 
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Figure 1(a) Figure 1(b) 
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Figure 2 


The DIS cross section with unpolarized beam and target can be written in the one 
photon exchange approximation as [ Q 


_ q-Mott t^E2{xb,Q‘^) 1 -f eR{xB, Q^) 
dxsdQ'^ EE' exs ^ + R{xb,Q‘^) ’ 


( 1 ) 


where e is a kinematic parameter and R is the ratio of longitudinal to transverse photon- 
hadron cross section, R = o^Iot = [(1 + ‘ix\rn?/Q"^) E 2 — 2xbEi\/[2xbEi]. The functions 
El and E 2 are called structure functions and dehned as 


Ei{xB,Q^) = ef'^W^,{xB,Q^) ( 2 ) 

with z = 1,2 and projection tensors ef' given in Ref. [ ||. The lF^y(xB,Q^) is the DIS 
hadronic tensor and is proportional to the square of the hadronic part in Fig. 1(a) or the 
imaginary part of the forward scattering amplitude in Fig. 1(b). The structure functions 
El and E 2 contain all complex hadronic interactions of the DIS cross section. 

The structure functions are locally dehned in DIS and can be directly extracted from 
measured DIS cross section via Eq. (0- Because DIS cross section depends on both 
hard scattering scale and soft momentum scale intrinsic to hadron wave function 
1/fm ~ Aqcd, structure functions Ei and E 2 are nonperturbative quantities. 

2.2. Parton distribution functions 

Parton distributions (or parton distribution functions), 0///j(x,/i^), are dehned as ma¬ 
trix elements of a pair of parton helds and are often interpreted as the probability densities 
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for finding a parton of flavor / within a hadron h of momentum p, with its momentum 
fraction between x and x + dx and it virtuality less than /i^ [ |^. For example, a quark 
distribution is given by 

= J . (3) 

Corresponding Feynman diagrams in momentum space are given by the type of cut-vertex 
diagrams in Fig. 2 [^. 

Parton distributions are universally dehned and in principle, in dependent of any spe- 
cihc physical process. Because of the hadron state, like the structure functions, parton 
distributions are nonperturbative quantities. However, unlike the structure functions, 
parton distributions are not direct physical observables. 
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2.3. Factorization 

Structure functions and parton distributions are related to each other via perturbative 
QCD factorization [P]. 

In QCD perturbation theory, the forward scattering amplitude in Fig. 1(b) can be 
represented in terms of an expansion of Feynman diagrams illustrated in Fig. 3. However, 
because of the soft momentum scale in the hadron state, Feynman diagrams in Fig. 3 
are not entirely calculable in perturbation theory. For example, the integration of loop 
momentum k for the leading order (LO) diagram can be expressed as 
-0.8cm 
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J,rkH{Q,k) {T^)T{k,p) 


which is completely dominated by the phase space near ~ (1/fm)^ ~ Aq^j^ because 
of the perturbative pinch singularity at ~ 0. Therefore, when -C Q^, the dominant 
contribution to the DIS cross section comes from the phase space where the active quark is 
“long-lived” relative to the time scale of hard collision, tc ^ 1/Q. It is such a “long-lived” 
parton state that separates long- from short-distrance effects in the cross section. 
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Figure 4 
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In terms of a power expansion of (A;^)/ Q^, the LO forward scattering amplitude in Fig. 3 
can be approximated as illustrated in Fig. 4, and corresponding contribution to structure 
function F 2 can be expressed as 

F2(XB,Q^}^XBY,eU,/l^(.XB) + o(^^^ . ( 4 ) 

If we neglect all high order diagrams and power corrections in 1/Q^, F 2 in Eq. (^) has the 
Bjorken scaling in and is the same as the prediction of the parton model. 



Figure 5 
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However, QCD is much richer in dynamics than Feynman’s parton model, and the rest 
diagrams in Fig. 3 will all contribute to the measured structure functions. Although there 
are inhnite diagrams, leading contributions to in a physical gauge come from the 
diagrams with a ladder structure, as illustrated in Fig. 5, and partons’ loop momentum 
integrations / d'^ki are dominated by the region of phase space near the perturbative pinch 
singularities, kf ~ 0. By neglecting the power corrections of {k‘f)/Q‘^, as illustrated in 
Fig. 5, the sum of all ladder diagrams in Fig. 5 can be factorized into a convolution of 
a short-distance part with all momentum scales of the order of and a long-distance 
nonperturbative parton distribution [ , 


F2(xB,Q^) = Y.f' 


/xb X 


X 


(()f/h{x,fi^) + O 


^QCD 


(5) 


where J2f runs over all parton falvors, the coefficient functions Cf are perturbatively 
calculable in a power series of as, and the parton distributions (j)f/h represent the leading 
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and universal part of the long-distance physics. All other nonperturbative contributions 
are powerly suppressed. 

If the scattering is hard enough for neglecting all power corrections, Eq. represents 
a fact that parton distributions, though theoretically dehned, can be experimentally mea¬ 
sured via the physical quantities like the structure functions. It also indicates that a test 
of QCD in high energy collisions requires a good knowledge of parton distributions. 


2.4. Global analyses of parton distributions 

Although parton distributions are nonperturbative, their dependence on the factoriza¬ 
tion scale /i^ is predicted by pQCD in a form of DGLAP evolution equations [ , 


d 




(h)) 


( 6 ) 


with calculable splitting functions Pi/j for /i^ 3> Aq,^^. To solve for the parton distribu¬ 
tions, we need a set of input parton distributions at /ig, which can only be extracted from 
experimental data. 

The global QCD analysis of parton distributions represents our effort to hnd a best set 
of universal parton distributions from all existing data. The analysis itself is an excellent 
test of QCD dynamics in the hard scattering, the factorization theorems, and the univer¬ 
sality of parton distributions. With the extracted parton distributions of a free nucleon, 
pQCD calculations with next-to-leading order (NLO) accuracy in are consistent with 
thousands of data points from more than dozen physical observables in hadonic collisions 
[ ^ Since the factorization theorems determine the absolute normalization for each 
observable, there is no need for any artihcial A'-factor. 


3. Nuclear Structure Fuuctious aud Partou Distributious 
3.1. Defiuitious 

The derivation of the DIS cross section in Eq. (|^) is independent of the details of the 
targets. Nuclear structure functions, F^^XByQ"^), extracted from DIS data on a nuclear 
target of atomic weight A, are dehned in the same way as that in Eq. (|^) with the 
hadronic tensor of a nuclear state. The Bjorken variable, xb = with an averaged 
nucleon momentum p = ^, and has a range from 1 to A. 

If the partons’ typical virtuality in a nucleus, Q^, and {kT)A -C xp, all deriva¬ 

tions in last section for a free nucleon state can be carried over for a nuclear state, and 
nuclear structure function, shares the same factorized relation, 

Fi(xs,Q^) = f + (7) 

j X X fJj j y j 

where the short-distance coefficient functions Cj should be the same as those for free 
nucleon and independent of A. In Eq. (0), nuclear parton distributions, are 

matrix elements of nuclear states with the same operators of nucleon parton distributions. 

If we can neglect the power corrections in Eq. (0), i.e., pQCD factorization theorems 
hold, we can extract nuclear parton distributions and their A-dependence from the mea¬ 
sured structure functions on nuclear targets. Because of the universality of parton dis- 
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tributions, the extracted A-dependence of nuclear parton distributions should also be 
universal and represent the internal properties of a nuclear wave function. 

-O.Tcm 



Figure 6 
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3.2. Nuclear shadowing 

Gluon distribution grows rapidly as momentum fraction x decreases. When xb ^ 0.1, 
a nucleus becomes a dense system of small-x gluons and the large number of soft gluons 
from different nucleons can all participate in the hard scattering, as illustrated in Fig. 6. 
Such coherent multi-gluon interactions suppress nuclear structure function F.f{xB,Q'^) 
in comparison to a sum of free nucleon structure functions. Such suppression is often 
referred as nuclear shadowing, and it increases as xb decreases and/or A increases. 



Figure 7 
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In order to understand the phenomenon of nuclear shadowing in a language close to 
pQCD factorization, let’s consider a leading twist (or leading power) gluonic contribution 
to the structure function F 2 , as illustrated in Fig. 7. If the active quark’s virtuality 
in a large nucleus, {k‘^)A ^ the multiple gluonic contributions are dominated by 
the phase space where the quark state is “long-lived”, which separates long- from short- 
distrance effects, and pQCD factorization holds. Unlike what shown in Fig. 5, nuclear 
quark distribution in Fig. 7 gets contributions from multi-gluon interactions from different 
nucleons. Since {k'^)A ^ such multi-gluon interactions are remote from the short 
distance hard scattering, and thus, are internal to nuclear parton distributions. 

If the factorization scale is within perturbative region, multi-gluon interactions in 
Fig. 7 lead to calculable corrections to DGLAP equations. The hrst corrections, in the 
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small-a: limit, are known and lead to a set of modified evolntion equations [ For gluon 
distribution in a nucleus, we have 

2 9 2 ^ (^sCa dx' , 2 ^ (asCA\^ f^dx' , , 2^^2,o^ 

f, = j I —l^GA^.^)] (8) 


9/i- 

where k is a known positive constant [0 


The second term slows down the rapid growth 
of gluon distribution at small-x due to gluon recombination. Although it is suppressed 
by a power of the effect of this nonlinear term to gluon distribution does not really 

suppressed by 1/u^ because of a simple fact that / 2 

Mo 


dfj,^ 


as /i 


00. 


m5 

According to pQCD factorization theorems, A-dependence of nuclear parton distri¬ 
butions are universal, and can come from two sources: (1) the nonlinear terms in the 
modihed evolution equations and (2) input nuclear parton distributions at Pq, which have 
to be extracted from experimental data. It was shown 10 that the nonlinear evolution 
in perturbative region cannot be a main source to generate the observed nuclear shadow¬ 
ing, and parton distributions exhibit signihcant shadowing at the initial scale of evolution 
However, for hxed small-a: the nonlinear evolution is sufficient to keep the shadowing 


ho 


going away very slowly as one increases /i^ 


even 


up to values as large as 100 GeV^. This 


is a solid prediction of perturbative QCD and has been confirmed experimentally. 

On the other hand, A-dependence of nuclear structure functions are not universal be¬ 
cause of the power corrections in Eq. (|^). 

3.3. Saturation — breakdown of pQCD factorization 

PQCD factorization for nuclear structure functions in Eq. (@) is only valid for a dilute 
nucleus. If the direct power corrections in Eq. (|^) and/or the power corrections to the 
modified evolution in Eq. (§) become important it is a signal that a dense system of 
partons is reached. When the power corrections are comparable to the leading power 
contributions, the conventional pQCD factorization breaks down. 

Although Eq. (|^) is really valid only when the nonlinear term is small compared to 
the usual evolution term. Nevertheless, Eq. (H) is an interesting equation in that the 
nonlinear term stabilizes the growth of normal evolution and leads to a limiting value for 
xGa^x^P^^) as a: —> 0. We can estimate roughly where this saturation sets in by finding 
where the gluon distribution loses its a:-dependence, and find xG{x,Q‘l) oc I^Ql/as{Qs) 
at a saturation scale More recent developments in understanding saturation can be 
found in Ref. 
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Fi{XB.Q^) 
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When xb —>• 0, all order gluonic contributions to the structure functions, as illustrated 


















in Fig. 6, are equally important. If the intrinsic parton virtuality in such a dense system, 
{k‘^)A ^ Qs ^ ^QCD) active partons are “short-lived”, and multiple gluonic contribu¬ 
tions to the structure functions in Fig. 6 or 7 might be calculated without introducing a 
“sea” quark distribution. The calculation can be carried out in the target rest frame as 


illustrated in Fig. 8, and the saturated nuclear structure function is given by [0 


Fq oc / Fbdz 


' 07 *—Q ) ^qq—A{b,Z,Qg') 


(9) 


where is the wavefunction for a virtual photon to go into a quark-antiquark pair 

of longitudinal momentum fractions 2 ; and 1 — 2 ;, 6 is a Fourier conjugate of relative 
momentum between q and g, and (Tqq-A reprensents a hadronic cross section between the 
qq pair and the nucleus. 

In this saturation regime, the pQCD factorization formula in Eq. (^) is not valid. 
What calculated in this saturation limit are nuclear structure functions including all 
power corrections, not the parton distributions, which, by dehnition, are twist-2 matrix 
elements and universal. 


3.4. Factorization in heavy ion collisions 

For a hard probe of scale Q in relativistic heavy ion collisions, there could be three types 
of multi-gluon interactions: (1) within individual ion, (2) between the hard scattering and 
one or both ions, and (3) between two ions, as illustrated in Fig. 9. 
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+ ... 


If the parton virtuality -C the hrst type soft gluon interactions are remote from 
the hard probe and internal to individual ion, and therefore, should have been included 
in nuclear parton distributions and do not interfere with pQCD factorization. The second 
type of multi-gluon interactions with physical polarizations are suppressed by powers of 
{k‘^)/Q‘^ because they link effects at two different scales. If soft gluon interactions with 
the hard part involve only one ion, a genralized factorization for leading power corrections 
should be valid for calculating this type of interactions 
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The third type of multi-gluon interactions is the most difficult one to remove. Such 
soft gluon interactions exchange informations between two ions, and therefore, have a 
potential to alter the parton distributions before the hard scattering takes place. Without 
the universality of parton distributions, pQCD calculations do not have real predictive 
power. Due to the gauge invariance of QCD and unitarity, contributions of the third type 
soft gluon interactions to a physical observable are suppressed by a power of {1/Q^) [ ^ 
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Therefore, pQCD factorization formalism for hadronic collisions can apply to hard 
probes in heavy ion collisions so long as the power corrections are relatively small, 

(PalA{Xa, ( 2 ) (j)b/B{Xb, 0 Hab + O (^ 0 ) 


where (8) repesents convolution over parton momentum fractions Xa and Xf,. One can also 
include nuclear thickness functions into Eq. (|^) to take care of the size effect of the ions. 
At this leading power level, A-dependence of an observable is completely determined by 
nuclear parton distributions. 

When the hard scale Q S, momentum fractions of the active partons can be very 
small. Soft gluon interactions between two ions, in addition to the interactions within 
individual ion, can lead to a much larger in heavy ion collisions than that in DIS. 
When Q is not too much larger than Qg, pQCD factorization fails, and normal concept 
of parton distributions cannot be applied to such a dense system. 

If the power corrections are signihcant, but not too large, pQCD calculations might 
be systematically carried out at the hrst power corrections. The predictions are often 
expressed in terms of multi-parton correlation functions 
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4. Global Analysis of Nuclear Parton Distribution Functions 
4.1. Methodology 

Predictive power of pQCD calculations of hard probes in heavy ion collisions relies on 
a set of universal nuclear parton distributions. Therefore, a reliable set of nuclear parton 
distributions is crucial for understanding RHIC physics. Any new physics and phenomena 
should represent the observations beyond what predicted by pQCD factorization. 

PQCD predicts the scale dependence of parton distributions in terms of evolution equa¬ 
tions. Before performing a global QCD analysis of nuclear parton distributions, we are 
required to choose: (1) evolution equations, (2) parameterizations of input nuclear parton 
distributions at a scale /ig, and (3) data sets with rich nuclear information. Then, we need 
to (1) evolve the input distributions to any other values of (2) calculate the theoreti¬ 
cal predictions by using the evolved distributions, and (3) compare the predictions with 
real data and calculate the The best set of nuclear parton distributions should give a 
minimum or a best £t to all existing data. 

At the leading twist level, A-dependence of all physical observables are completely 
determined by the nuclear dependence of the parton distributions. Deriving a reliable 
A-dependence of parton distributions is a most crucial part of the global analysis. There 
are two possible approaches to derive the A-dependence: (A) DGLAP evolution with 
all A-dependence included in the input distributions at pg and (B) modihed evolution 
equations of the type in Eq. (||) with nonperturbative A-dependence included in the 
input distributions at /ig and perturbative A-dependence generated from the evolution. 


4.2. Exsiting work 

Two groups have been doing global analysis of nuclear parton distributions. Eskola, 
Kolhinen, Ruuskanen, and Salgado (EKRS) produced EKS98 package of nuclear parton 
distributions [ Hirai, Kumano, and Miyama (HKM) derived several sets of nuclear 


parton distributions from extensive DIS data [ . 
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Both groups used LO DGLAP evolution equations and adopted the approach (A) for 
including nuclear dependence. They dehne the input parton distributions as 


0/M( 2 :, ho) = ho) ho) 


( 11 ) 


with are known free nucleon parton distributions. EKRS used both CTEQ- 

LO and GRV-LO free nucleon parton distributions for producing EKS98, while HKM 
used MRST-LO parton distribution for its nuclear parton distributions. Two groups used 
different parameterizations for Rf{x, /Xq). The goal of the global analysis is to extract the 
ratio Rf{x,fil) that represents a best £t to all data used in the analysis. 

For data selection, two groups used different data sets to produce their published nuclear 
parton distributions [11^, [T^. EKRS used both DIS and Drell-Yan data on nuclear targets. 


In particular, Drell-Yan data from Fermilab E772 experiment help to £x relative nuclear 
effects in valence and sea and make the ratio for sea quarks, Rs{x,fiQ) < 1 at medium x. 
In addition, EKRS used NMG data on Q^-dependence of Tf “/ F 2 that provide constraints 
on R^{x,nl). On the other hand, HKM used only DIS data. 

Detailed comparison between data and theoretical calculations using EKRS nuclear 
parton distributions and other available parameterizations can be found in Ref. [01- 
Similar comparison for HKM distributions can be found in Refs. 
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Both published 

EKRS and HKM nuclear parton distributions can £t their selected data sets very well. I 
summarize a few key differences here: (1) there are large differences in input distributions 
at /ig = 2.25 GeV^; (2) the default set of nuclear parton distributions in HIJING seems to 
give a too strong A-dependence; and (3) HKM distributions appear to predict a different 
Q^-dependence for NMG data on the ratio F 2 and a different shape in X 2 dependence 
for E772 Drell-Yan data. 

Notice that both of these data were not included in the HKM’s original global analysis. 
Since DIS data are only sensitive to the sum of quark and antiquark distributions, it should 
not be surpprised that the separation of sea and valence quarks in HKM distributions is 
not very well constrained. Because LO structure function F 2 does not have an explicit 
dependence on gluon distribution, gluons could not be well constrained by DIS data alone. 
These two ambiguities are reduced in EKRS approach because of the usage of Drell-Yan 
data and NMG data on Q^-dependence of F^'^/F^. 


5. Summary and Outlook 

Much of the predictive power of pQGD calculations relies on the factorization and uni¬ 
versality of parton distributions. PQGD calculations with next-to-leading order accuracy 
in as have been consistent with almost all existing data from hadronic collisions without 
a need for an artihcal iP-factor. 

PQGD factorization theorems can be applied to hard scattering involving nuclei, if the 
partonic system covered by the interaction volume is a dilute system, in which partons 
are “free” or “long-lived” at the time scale of hard collision. If the leading twist pQGD 
factorization theorems hold, A-dependence of all observables are completely determined 
by the universal A-dependence of nuclear parton distributions. Any signihcant deviation 
from the predicted A-dependence is a signal of breakdown of the conventional pQGD 
factorization and new physical phenomena. Therefore, precise A-dependence of nuclear 
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parton distributions are very important for understanding hard scattering signals in heavy 
ion collisions at RHIC and the LHC energies. 

Nuclear dependence of the input distributions is a main source of the H-dependence of 
nuclear parton distributions and has to be extracted from better data. At high energy 
or small-x, perturbative A-dependence via the modihed evolution equations might be 
important too. Two sets of nuclear parton distributions from LO global QCD analysis 
are available. In order to reach the same level of precision tests that achieved in collisions 
with free hadrons, we need NLO nuclear parton distributions with better precisions on 
the A-dependence. 
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